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Abstract: The Tibet air siiower array, which has an effective area of 37,000 square meters and is located 
at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few 
TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to 
cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic- 
ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in 
our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte 
Carlo simulation. 
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Introduction 



Although SNRs are theoretically considered to be 
the most plausible candidates for acceleration of 
cosmic -ray hadrons up to PeV energies, no ob- 
servations have succeeded in identifying them so 
far. Since accelerated electrons have difficulty pro- 
ducing very high-energy gamma rays with ener- 
gies above 100 TeV via bremsstrahlung or inverse 
Compton scattering, it can be an effective way of 
obtaining clear evidence for hadronic acceleration 
to detect high-energy gamma rays above 100 TeV 
generating via the decay of neutral pions produced 
in interactions of accelerated hadrons with ambient 
material, e.g. molecular clouds. In other words, 
study of the energy-dependent morphology of the 
sources in the 100 TeV region (10-1000 TeV) is a 
clue to disentangle their acceleration mechanism. 

The HESS group reported in 2005 on the discovery 
of 14 new TeV gamma-ray sources in the southern 
hemisphere and many of them have a harder energy 
spectrum (indices: -1.8 to -2.8) [1], which implies 
that hadrons should be accelerated there. No defi- 
nite claims have been made yet, however, because 
it is difficult for HESS to measure the sources 
above ~10 TeV. Also, many of the 14 sources are 
faint in other wavelengths. This strongly suggests 
that an unbiased survey irrelevant to observations 
in other wavelengths using a wide field-of-view ap- 
paratus, e.g. an air shower array, is needed. 

The Tibet Air Shower Experiment 



30 m 



Figure 1: Schematic view of the Tibet AS+MD ar- 
ray. Open squares and open circles represent the 
surface scintillation detectors that compose the Ti- 
bet AS array. Note that the AS array drawn here 
is upgraded from the current version so that its ef- 
fective area becomes 50,000 m^ by modifying the 
configuration of the scintillation detectors. Filled 
squares show the proposed Tibet MD array 2.5 m 
underground. 

by the Milagro group [6J. We can not draw any 
clear conclusion on the region, however, because 
the AS array at present is unable to discriminate 
gamma rays from hadrons in the multi-Te V energy 
region. 



The Tibet air shower experiment has been success- 
fully operating at Yangbajing (90°31' E, 30°06' N; 
4,300 m a.s.l.) in Tibet, China since 1990. Being 
expanded several times, the Tibet air shower (AS) 
array now consists of 789 plastic scintillation coun- 
ters placed on a 7.5 m square grid with an effective 
area of 37,000 m^, and is detecting high-energy (> 
a few TeV) cosmic gamma rays and hadrons. 
Using the Tibet AS array, we so far have observed 
a new cosmic-ray anisotropy in the Cygnus region 
at multi-TeV energies [2] as well as TeV gamma 
rays from the Crab [3], Mrk501 [4] and Mrk421 
[5]. Because the anisotropy that we observed in the 
Cygnus region is relatively narrow, it may favor the 
existence of a gamma-ray component as claimed 



The Tibet Muon Detector Array 

The currently proposed configuration of the Tibet 
MD array is shown in Figure 1 . It is composed of 
12 pools, each of which consists of 16 cells. Each 
cell is a waterproof concrete tank which is 7.2 m 
wide X 7.2 m long x 1.5 m deep in size. Two 
20 inch-in-diameter photomultiplier tubes (PMTs, 
Hamamatsu R3600) are put on its ceiUng, facing 
downwards. Its inside is painted with white epoxy 
resin to waterproof and to efficiently reflect water 
Cherenkov light, which is then collected with the 
PMTs. The MD array is set up 2.5 m underground 
(~19 radiation lengths) in order to detect the pene- 
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trating muon component of air showers, suppress- 
ing the electromagnetic one. Its total effective area 
amounts to 9,950 for muon detection with the 
energy threshold of approximately 1 GeV. 

Results and Discussions 

We performed detailed Monte Carlo simulations 
regarding air showers and the responses of soil, 
the AS and MD array, which are described in [7]. 
Based on the number of collected photoelectrons 
by the PMTs of the MD array, background hadrons 
can be rejected by selecting muon poor events. 
Around 100 TeV, the number of hadron-induced 
events are suppressed down to 0.01 % or less, while 
gamma-induced events are retained by more than 
83%. 

Table 1 summarizes the comparison of perfor- 
mance between HESS and the Tibet ASh-MD ar- 
ray. One can see the advantages of an air shower 
array, i.e. a wide field-of-view (—1.5 sr) and high 
duty cycle (—90%). Although the angular and en- 
ergy resolution of the ASh-MD array around 100 
TeV are two times worse than that of HESS, the 
difference is not so critical. Taking RX J 17 13 as 
an example, a famous source in the southern sky 
with radius of 0.6° and spectral index of — —3, the 
ASh-MD array could extract meaningful physical 
quantities. 

The sensitivity of the Tibet AS+MD array to a 
gamma-ray point source deduced from the simula- 
tion is shown in Figure 2. Its 5tT sensitivity in one 
calendar year will reach 7% and —20% Crab above 
20 and 100 TeV respectively, and surpass the exist- 
ing lACTs above 20 TeV. Furthermore, it may sur- 
pass the sensitivity of the next generation lACTs 
above 40 TeV. Since no extensive search has been 
conducted with an apparatus with sensitivity com- 
parable to HESS in the northern sky, the AS+MD 
array could discover as many (—10) unknown TeV 
gamma-ray sources as HESS did in [1]. 

As shown in Figure 3, the Tibet AS+MD array will 
have more advantage in observing spatially diffuse 
gamma-ray sources, e.g. MGRO J2019+37, dis- 
covered by the Milagro experiment [8]. 

In brief summary, among the known sources in 
our field-of-view. Crab, TeV J2032+4130, HESS 
J1837-069, Mrk 501, Mrk 421, MGRO J2019+37, 



Table 1: Comparison of performance between 
HESS and the Tibet AS+MD aiTay 
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Figure 2: The attainable integral flux sensitivity 
(5(7 in one calendar year) of the Tibet AS+MD ar- 
ray to a point-like gamma-ray source, together with 
the sensitivity of HESS and some other future ex- 
perimental plans. 
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MGRO J2031+41 and MGRO J 1908+06 [6] are 
sufficiently detectable. HESS J1834-087, LS 1+61 
303, Cas A, and M87 are marginal. 



MGRO J201 9+37 (sigma = 0.32 deg., 1 .7 radius analysis) 
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Figure 3: 5a sensitivity of the Tibet AS+MD ar- 
ray in one calendar year (red) to a diffuse gamma- 
ray source [8]. The blue line is expected from the 
Cygnus region in the galactic plane. 
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